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ABSTRACT 



X 

■ We report here the first study of proper motions of fast filaments in the 

young, oxygen-rich supernova remnant G292.0+1.8, carried out using a series of 
[O III] 5007 A emission-line images taken over a period of more than 21 years. 
Images taken at seven epochs from 1986 to 2008, all from telescopes at the 
Cerro Tololo Inter- American Observatory, show oxygen-emitting filaments, pre- 
sumably ejecta fragments, throughout most of the remnant. We have measured 
the proper motions for 67 discrete filaments through two-dimensional correlations 
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between images from different epochs. While the motions are small, mostly 20 
to 100 masyr"^, they are nevertheless measurable through a robust technique of 
averaging measurements from many epoch pairs. The data are qualitatively con- 
sistent with a frcc-cxpansion model, and clearly show systematic motions outward 
from a point near the center of the radio/X-ray shell. Global fits using this model 
indicate an expansion center at R.A. (2000.0) = ll''24™34.4^ Dec. (2000.0) = 
— 59°15'51", and a kinematic age of 2990 ± 60 years. The young pulsar PSR 
J1124-5916 is located 46" southeast of the expansion center. Assuming that it 
was launched by the supernova, we expect the pulsar to be moving southeast- 
ward at 16 masyr"^, or a transverse velocity of 440 kms~^. We find the fastest 
ejecta along an axis oriented roughly N-S in the plane of the sky, suggesting that 
a bipolar explosion produced G292.0-I-1.8, as appears to have been the case for 
Gas A. 

Subject headings: ISM: individual (SNR G292.0+1.8) — ISM: kinematics and 
dynamics — supernova remnants 

1. Introduction 

Massive stars produce copious quantities of oxygen and other low-Z elements during 
their hydrostatic evolution, and these arc ejected when the stars end their lives in spectacu- 
lar fashion as core-collapse supernovae (SNe). After the pyrotechnics of the explosion fade, 
the ejecta coast outward at velocities of up to several thousand kms~^, and the fastest of 
them interact with the surrounding circumstellar and/or interstellar medium and produce an 
outward propagating blast wave. The interaction also results in the dramatic slowing of the 
ejecta and in a reverse shock that propagates inward, gradually encountering slower-moving 
ejecta and exciting them, until eventually the shock dissipates as it ncars the center. As 
dense knots of ejecta encounter the reverse shock, they are excited and become visible, ra- 
diating emission lines characteristic of their composition — oxygen especially since this is the 
most abundant element in the outer core of typical progenitors. The resulting oxygen-rich 
filaments — with a spectrum dominated by oxygen lines, showing little or no hydrogen, and 
with typical velocities > 1000 kms~^ — may be seen for a few hundred to a few thousand 
years. Individual filaments last for a much shorter time, as their material is evaporated 
through the interaction with the reverse shock to enrich the hot. X-ray-emitting plasma 
that fills most of the shells of these young supernova remnants (SNRs), producing an X-ray 
spectrum with strong lines from elements including Ne, Mg, Si, S, and others in addition to 
O. Another product of a core-collapse SN is a compact remnant — a neutron star or perhaps 
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a black hole — that may also be visible as a point, possibly pulsing, X-ray source, sometimes 
surrounded by a synchrotron-emitting pulsar-wind nebula. Observationally, objects display- 
ing these fast, 0-rich optical filaments have become known as oxygen- rich SNRs, a small 
class including only Cas A, Puppis A, and G292.0-I-1.8 in the Galaxy plus 4 to 6 others in 
the Magellanic Clouds and other nearby galaxies. 

Of these, the only one to show all the above characteristics expected from a core-collapse 
supernova is G292.0+1.8 (hereafter simply G29 2 ). G2 92 was first discovered and denoted 
as MSHll-5^ in the radio survey bv lMills et al.l (119611 ), and on the bas i s of it s non-thermal 
spectrum was identified as a SNR by iMilnel Jl969h and IShaver &: GossI (119701) . Much more 
recently, the most detailed radio images of G292 are by iGaensler fc Wallace! (120031 ) . who 
found that the SNR has a diameter of 8', and who also give what is probably the most 
reliable distance estimate: 6.2 ± 0.9 kpc, based primarily on the H I absorption profile. 
Throughout this paper we will scale all distance-dependent quantities to = d/{6 kpc). 



X-rays from this source were first identified in HEAO-1 data by Share et al. 
and its properties were investigated from the Einstein O bservatory by 
who found strong X-ray lines from Mg, Si, and S, and by iTuohy et al.l ( 



Clark et al. 



1978 ) 



1980) 



1982h, who noted a 



prominent central ring, possibly suggesting a rotating Type II SN precursor. Mo st recently, 



studies of G 292 from Chandra have produced spectacular high-resolution images (jPark et al. 



2OO21 . 120071 ) . The X-ray spectrum is dominated by K-shell lines of O, Ne, Mg, Si, and S. Non- 
equilibrium ionization analyses require large enhancements (relative to solar) in abundances 
for O, Ne, and Mg, with lesser enhancements for Si, S, and Fe (iHughes fc Singhl Il994l : 



Gonzalez &: Safi-Harbl 120031 ) . By comparing the inferred abundances with the integrated 
yields predicted by models for core-collapse supernovae, these groups estimated a progenitor 
mass in th e rang e ~ 25 — 40 Mq. In a more detailed study of individual X-ray features. 
Park et al.l (120041 ) found that different knots have very different compositions and suggested 
that these represe nt clumps of eject a from different zones of the progenitor, and in a much 
deeper ACIS image iPark et al.l (120071 ) have found evidence for asymmetries in the distribution 
of (at least) the oxygen ejecta. These latest images show that the X-ray shell extends to the 
south, with a full extent of 9.6' (N-S) x 8.4' (E-W), or 16.7 x 14.7 4 PC 

The original Chandra image of G292 also reve aled a compact cent ral source, surrounded 
by what app e ared to be a pulsar wind nebula (IHughes et al.ll200ll ). Shortly thereafter. 



Camilo et al 



Hughes et al. 



(I2OO2 



discovered the radio pulsar PSR J1124-5916 within or very near G292. 
(l2003l ) then showed that the compact X-ray source is pulsed with the same 
period, confirming that it must be the compact remnant of the star that produced G292. 
The period of 135 ms and spi n-down age of 29 00 yr are roughly consistent with the age of 



the G292 SNR, estimated by IChevalierl (120051 ) as 2700-3700 yr based on properties of the 
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pulsar-wind nebula. 



Goss et al.l (Il979l ) discovered optical emission from G292, consisting solely of a group 



of filaments near the center of the radio shell with a spectrum sho wing strong lines from O 
and Ne, and no lines from H, N, or S that are more typical of SNRs. iMurdin fc ClarkI (119791 ) 
showed that the brighter filaments have high radial velocities, —650 < Vrad ^ +1380 kms~^, 
suggesting that they are undecelerated (or minimally decelerated) ejecta from the supernova 
explosi on. Together the s e resu lts provided the first identification of G292 as an oxygen-rich 
SNR. iMurdin &: ClarkI (119791 ) estimated the remnant's age as r < D j lS.v ~ 1800 yr, 
based simply on a scaling argument using a diameter D '2! 3.5 pc, the extent of 
the central radio peak where they had observed the optical filaments, and their m easured 
radial velocity range of 2030 kms~^. More limited spectra bv Ivan den Berghl (19791) showed 



a velocity dispersion only slightly smaller, and further spectra by iBraun et al 
a similar kinematic estimate for the age, r ~ 2000 yr,^ 



(Il983l ) led to 



Ghavamian et al.l (120051 ) . in an extensive kinematic study using the Rutgers Fabry-Perot 



imaging spectrometer, reported [O III] A 5007 emission features from the central 7.5' of G292 
throughout a velocity range of at least —1440 < v^ad ^ 1700 kms~^ (the full range covered in 
their etalon settings). They found the distribution of 0-rich knots to be quite asymmetric, 
with virtually all of the knots north of the (radio-determined) shell center being blue-shifted, 
while toward the south ther e are both red- and blu e-shifted knots. From a plot of radial 
velocity vs projected radius, iGhavamian et al. J2005h found fair agreement with a model in 
which the knots are distributed around an expanding shell. Fitting this distribution with 
a velocity-radius ellipsoid, they found an ejecta shell velocity Vej ~ 1700 kms~^, centered 
about a systemic velocity of Vc ~ +100 kms~^. Assuming a spherical shell geometry led to 
an estimated age of r = R^j/Vej = (3000 — 3400)^6 yr, where the range refiects both the 
width of the velocity-radius ellipsoid and models us ing somewhat different assu mptions. The 



greate r age, by almost a factor of 2 compared with iMurdin fc ClarkI (119791 ) and iBraun et al. 



(119831 ). is now quite consistent with the pulsar spin-down age and the estimate based on the 
PWN pr operties. The change iii the optically-based estimate is due almost entirely to the 
fact that IGhavamian et al.l (120051 ) observed (fainter) 0-rich filaments out to a radius roughly 
twice that found in the earlier studies. 



The full optical extent of G292 was shown for the first time by IWinkler fc Long] (12006 



paper 1), who found 0-rich optical knots throughout most of the 8' shell seen in radio and 



^Both Murdin fc Clark ( 1979t) and Braun et al. (1983) used a distance of 5.4 kpc and hence found slightly 
smaller ages. We have simply scaled their results to a uniform distance of 6 kpc. Since they detected only 
the brightest optical filaments, near the center of the shell, their estimates for the diameter and scaled age 
are much smaller than more recent ones. 
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X-ray images, with a few knots extending even beyond th e the primary sh ell and into the 



southern extension seen in the recent deep Chandra image (iPark et al.l 120071 ). In paper 1 we 
also found a few faint knots showing [S II] AA 6716, 6731 emission in addition to the oxygen 
lines, the first optical evidence for the presence of 0-burning products in the ejecta. Most 
of the outer knots exhibit radially-oriented pencil-like morphologies, suggesting an origin as 
Rayleigh- Taylor fingers. 

In this paper we present our study of the proper motions of oxygen-rich knots in G292, 
the result of a series of CCD images from 1986 to 2008. These clearly show the systematic 
motion of the knots outward from a point near the center of the outer radio shell. Global fits 
of a free-expansion model lead to what we believe to be the most reliable age measurement 
for the SNR, and indicate that the center is 46" from the pulsar PSR J1124-5916. We discuss 
briefly the three-dimensional distribution of ejecta, which indicates a clear asymmetry and 
suggests that the progenitor explosion may have been a bipolar one. 



2. Observations and Data Reduction 

Emission-line images used in this study were taken in seven separate observing runs 
on various telescopes at the Cerro Tololo Inter-American Observatory (CTIO) from 1986 
through 2008. A journal of the observations is given in Table 1. The runs from 1986 and 
1991 used CCD chips with a small field of view that covered only the central region of G292, 
as noted in the table. Those from 1999 onward all used the 0.9m telescope and Tektronix 
2K#3 chip, and covered a field 13.7' square, easily encompassing the entire remnant, at 
a scale of 0.401" pixel"^. All the observations used narrow-band interference filters, with 
characteristics also listed in Table 1. Note that while some filter nominally designed for 
[O III] A 5007 was used in each case, we in fact used three different such filters, and even 
for the same filter the effective bandpass is different when used in beams converging with 
different focal ratios. In some cases, especially the runs on the 0.9m in 1999 and 2000, the 
filter was one designed for a beam faster than the 0.9m's of f/13.5; as a result some of the 
most blueshifted filaments are barely visible in the 1999 and 2000 images. In addition to 
images in [O III] A 5007, we also obtained images in a continuum band at each run from 
1991 onward, so we could subtract away most of the myriad stars that pervade the crowded 
Galactic field of G292 and make it easier to measure the motions of the filaments. 

On the 1986 run we obtained only a pair of [O III] images, but on all the subsequent 
runs we obtained three or more individual frames in both line and continuum filters, dithered 
by several arcsec between frames in order to remove bad columns and other systematic 
effects on combining the images. All images from all epochs were processed using standard 
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IRAF@ procedures for bias-subtraction and flat-fielding, the latter based on well-exposed 
dome or twilight sky flats. During each set of obs ervations we obtaine d several images of 



spectrophotometric standard stars from the list of iHamuy et al.l (119921 ) . and used these to 
flux-calibrate our images. 

In order to place all the images on a common wo rld coordinate systern (WCS), we used 



about 300 astrometric stars from the UCAC2 catalog (jZacharias et al.ll2003l ) (although many 
of these stars fell outside the fields of the earlier data). The stellar positions were corrected 
for proper motions to the appropriate epoch using the software provided with the UCAC2 
catalog. We then calculated a precise WCS for each frame using the IRAF tasks ccfind and 
ccmap to find the centroid of the reference stars and then fit them to a common projection. 
In the fitting process, we allowed for distortions from a standard tangent-plane projection of 
up to third order (using the "tnx" projection in IRAF). With the 0.9m frames we typically 
fit the positions of over 200 stars, with an rms dispersion of < 50 mas in both directions. 
For the earlier frames there were fewer reference stars, but equally good fits. 

We then defined (arbitrarily) a standard coordinate system: a tangent-plane projection 
centered at R.A. (2000.) = ll'^24™31.0^ Decl. (2000.) = -59°15'30.0", with a scale of exactly 
0.200" pixel"^, and each individual frame was transformed to the standard one using the 
IRAF task wregister. All the transformations used bilinear interpolation and double-precision 
arithmetic. The new scale was chosen as significantly finer than the pixel size for any of our 
images, while maintaining images of manageable size. At this point all the individual images 
from each run were combined to obtain a single image at each wavelength for each epoch. 

Precise quantitative measurement of the motions of knots is simplified greatly if we 
remove the stars, so that an individual knot can be isolated for the correlation analysis de- 
scribed in the next section. For this we used the matched continuum images, also transformed 
to our standard system and combined in exactly the same way as the [O III] imagesEI The 
continuum image for each epoch was appropriately scaled and subtracted from the emission- 
line image, judging by eye what scaling factor did the best job of removing the stars. If 
the stellar profiles of the emission-line and continuum image were not a close match, we 
used the IRAF task psfmatch, which determines a kernel to convolve with the sharper of 
the two images in order to match the profiles for stars in the less-sharp image. Happily, the 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by AURA, Inc., 
under cooperative agreement with the National Science Foundation 

■^For the 1986 image, where we had no matched continuum, we instead used a combination of continuum 
frames from later runs, suitably scaled and psf-matched. While less than ideal, this nevertheless did a good 
job at removing most of the faint stars. 
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continuum images had slightly narrower profiles in most cases, so we usually did not have 
to blur the resolution of the images showing the knots themselves. 

Once all the images were precisely registered on a common coordinate system, proper 
motions from one epoch to another can be seen, either from blinking images from different 
epochs, or from taking the difference between the two. In Fig. 1 we show an [O III] image of 
the entire G292 remnant, and in Fig. 2 we show enlarged versions of the difference images 
between different epochs for the locations indexed on Fig. 1. From such visual inspection of 
the images, we selected 67 individual knots or filaments for analysis. In order to be selected, 
filaments had to be: (1) sufficiently isolated that a rectangular box could be placed around 
the filament, enclosing virtually all of the emission and some surrounding background sky 
(this was necessary for our correlation analysis, §3.1); and (2) have no significant trace of 
residual stars on the filament itself in the continuum-subtracted images. (Residual stars 
in the surrounding background sky could be edited out manually, but we were afraid that 
attempting this on the filaments themselves could bias the data.) Some filaments were 
outside the field of view in images from epochs 1986 and/or 1999, or blueshifted out of 
the pass band in the 1999 and 2000 images, but all appear in images from at least three 
epochs: 2002, 2006, and 2008. We found no filaments that either appeared or disappeared 
ov er the 21-year span of our observ ations, unlike the situation with the much younger Gas 



A (IKamper &: van den Berghl 119761 ). 



3. Proper-Motion Measurements 

All the filaments used in our analysis were selected explicitly because of their strong 
[O III] A 5007 emission and absence of Ha, facts that strongly suggest these are fragments 
from the progenitor core that have interacted little if any with the interstellar medium since 
being launched by the explosion. It is thus reasonable to expect that the fragments we now 
see as optical filaments have been coasting with near-constant velocities since the explosion, 
a model which we investigate in detail in §4. For the moment, however, we describe in this 
section how we obtain the best values for the motions of the individual filaments. The only 
assumption is that the motion of each filament has been at essentially constant velocity over 
the 21-year history of our observations, a period that is less than one percent of the age of 
the SNR. 

The most straightforward approach to the problem might be to measure the position 
of each filament directly as a function of time using each of the epochs whe re we have 



observations. This is essentially the approach taken by iThorstensen et al.l (120011 ) for Gas A 



who then did a linear fit to measure the proper motions (see §3.3 for more on this approach). 
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For G292 this approach has two difficulties: most importantly, virtually all the knots and 
filaments have an irregular morphology, making it difficult to consistently define a precise 
center. This problem is compounded by the second one; namely, that individual filaments 
may change somewhat over time. (For exampl e fresh material raay enc ounter the reverse 



shock and light up, while other material fades.) iThorstensen et al.l (120011 ) encountered both 
these problems in Gas A, but the impact there was less important since the total baseline 
was 48 years, almost 20% of the age of the remnant, compared with less than 1% in the case 
of G292. 

In our case the motions are quite subtle; no filament shows a total displacement > 2" 
over the course of our observations, and typical epoch-to-epoch displacements are much 
less — certainly much less than the precision to which we could define the absolute position 
of a filament at any epoch. (The cases shown in Fig. 2 are among the most extreme.) But 
it is not necessary to precisely measure a filament's position on an absolute frame; precisely 
measuring its shift in position from one epoch to another, combined with an approximate 
measurement of its absolute position at some reference epoch, is sufficient for our purposes. 
For the relative position measurements we used a two-dimensional correlation technique 
described in the next section. 



3.1. Individual Measurements: Two-Dimensional Correlations 

To measure the displacement of an individual filament a from epoch i to epoch j, we 
first clipped out identical small rectangular sections from the aligned, continuum-subtracted 
images at the two epochs. Typical sections measured 30 to 100 pixels (6" to 20") in both x (E- 
W) and y (N-S) directions (exactly the small rectangles shown in Fig. 1), and were chosen 
to isolate an individual filament. We then used the IRAF task xregister to calculate the 
displacement, {Axa,i,j, Aya,i,j) from epoch i to j that gave the best match for the filamentary 
emission at the two epochs. The xregister task is a two-d i mensio nal implementation of 



the cross-correlation technique described by iTonry &: David (119791 ). Since the shifts are 



all < 2" (10 pixels), we used a discrete search over a small region, followed by a centroid 
fit to determine the best fractional-pixel shift in both dimensions. The cross-correlation 
technique does not give a reliable uncertainty estimate for an individual measurement; we 
use the scatter among measurements derived from multiple epoch pairs to obtain an overall 
uncertainty for each filament, as described in the next section. 
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3.2. Multiple-epoch Measurements: Averages, Weighting, and Uncertainty 

Our next step was to combine all the measurements over multiple baselines to give a 
mean value for the proper motion {fJ^xa^ fJ'ya) each of the 67 filaments, along with an 
uncertainty estimate. We did this independent of any model assumptions about the long- 
term kinematics of the filaments; i.e., we did not assume a free-expansion (or any other) 
model. In this process we analyzed each transverse dimension, x and y, independently; in 
what follows we will describe the analysis in the x direction; that for the y is, of course, 
identical. 

For each filament a, we calculate the mean proper motion //^^ by taking an average of 
the measurements of Axa,i,j, weighted by the square of the baseline Atij = tj — tf. 

V ^^a,»,j /\f2 VAr ■ At- ■ 
_— _ i,j — M. (^\ 

1,3 1,3 

where the shorthand notation ^ is used to indicate a sum taken over all the unique baselines, 

i,3 

ZlZl- (To see that At"^ is the correct weighting, consider that in combining data values 

j i<3 

with difi^erent known signal-to-noise, one weights as l/o"^, or [S/NY- For our correlation 
measurements, the noise is approximately the same for any epoch pair, but the signal, i.e., 
the displacement, increases linearly with the baseline. Hence S/N gc Atij, and Atfj is the 
best weighting factor.) For filaments that appear in observations from all seven epochs, there 
are 7x6/2 = 21 distinct basehnes to be considered. For those filaments that are present in 
images from fewer epochs (because they he outside the small fields observed in 1986 and/or 
1991, or because they were blueshifted outside the filter bandpass in 1999 and 2000), there 
are correspondingly fewer baselines. We shall denote the total number of baselines included 
in the average for each filament as A^. 

The advantage of using this multiple-baseline approach, compared to a simpler approach 
we describe briefiy in §3.3, is that we have a distribution comprising far more measurements 
whose variance we can calculate in order to estimate the uncertainties. The (unbiased) 
variance S^^ in the sample of measurements i^xodj with weights At'^j is 

Si - ^ " ^^^'^ 

1,3 I 1,3 
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and the uncertainty cr{fixa) iii the mean measurement fix^ is best estimated by 



^(/^.J' = -T^ -2SI (3) 



e.g., iBevington fc Robinsonll2002l ). 



In Fig. 3 we show a plot of individual proper-motion measurements, ^ix^ij ^"^^ l^vaiji 
as a function of time difference Atj j for all baselines, for three prominent filaments (the 
same ones highlighted in Fig. 2). Table 2 gives the measured values of positions and mean 
motions ^xai f^ya (with uncertainties) for all 67 filaments. We measured the positions in 
the 2002 image (the one with the best seeing) and extrapolated back to epoch 2000.0. We 
present the same data graphically, as vectors representing the proper motions extrapolated 
1000 years into the future at the present rate, superimposed on an image of G292 in Fig. 4. 
This figure clearly shows, at least qualitatively and as projected onto the plane of the sky, 
that the filaments are moving outward from a point near the center of G292, and that the 
ones farthest from the center are moving the fastest. We explore such a free-expansion model 
quantitatively in §4. 

We note that the measured proper motions are small, only 20 to 133 mas yr~^, with only 
4 of the 67 (all in the extreme south) having motions > 100 masyr"^. This is far smaller 



than for the closer and much younger oxygen-rich SNR Gas A, where iThorstensen et al. 



(120011 ) measured filamentary proper motions of ~ 160 to 800 masyr ^. 



3.3. Alternative Linear-Fit Approach 

An alternative to the multiple-epoch averages just described would be simply to plot 
the (absolute or relative) position of a filament as a function of time, and determine its 
mo tion from the slope of a linear fit. This is essentially the approach used for Gas A 



by IThorstensen et al.l (120011 ) . who measured absolute filament positions on a standardized 
reference grid. We have tried the same technique using the relative positions, measured as 
described in §3.1. In Fig. 5 we show the results for the same three prominent filaments as 
those highlighted in Figs. 2 and 3. 

This linear-fit approach provides a perfectly reasonable measure for the proper motions, 
but it does not extract the maximum amount of information from the data at hand. For 
a filament with observations at all seven epochs, the linear-fit method uses only six epoch 
pairs, compared with 21 for the multiple-epoch averages described in §3.2. Not surprisingly. 
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the uncertainties with the hnear-fit method are generally larger than for the weighted multi- 
epoch averages (§3.2). Nevertheless, we went on to fit both versions of the data using the 
free-expansion model, and obtained similar results (§4). 



4. Results: Global Kinematics and Free-Expansion Model 

The simplest model for the overall kinematics of the filaments in G292 is that the 
filaments we see today have been coasting with constant velocity since being launched by 
the explosion, and only recently rendered visible through interaction with a shock. Such a 
model has been successfu lly applied to th e other two oxygen-r ich SNRs in the Galaxy, Gas A 



(IThorstensen et al.ll200ll ) and Puppis A (jWinkler et al.lll988l ). Here we will use it for G292. 



For undecelerated motion, a filament now located at position r = (x, y) is expected to 
have proper motion {iJ^xai l^ya) given by 

_^ dr (r-ro) {x^ - Xq) iVa - Vo) 

- ^ = — = = ' (4) 

where fq = (xq, yo) gives the coordinates of the center of expansion (presumably the location 
of the progenitor at the time it exploded), and r is the age. Of course {xo,yo) and r are 
parameters of the model. 

Considering the ensemble of all 67 filaments, we have plotted components of the proper 
motion {fi^, function of {x,y), with the results shown in Fig. 6. The error bars 

represent the uncertainties as calculated in §3.2 and given in Table 2. It is obvious that 
these data are generally consistent with constant-velocity expansion, so we have performed 
linear fits to obtain values for xq (the x-intercept) and r (the inverse slope). The results of 
this analysis are: 

R.A. (2000.0) = ll''24™34.4^ Dec. (2000.0) = -59°15'51", r = 2990 yr. (5) 



Estimating the uncertainties to attach to these values requires more than a formal error- 
propagation analysis. Evaluation of for the linear fit (Eq. 4) to the data of Table 2 (and 
Fig. 6), we find an unacceptably large value: 5.3 per degree of freedom. (There are 67 x 
2 coordinates — 3 parameters = 131 degrees of freedom.) This should not come as a surprise, 
and does not mean we should reject the model entirely. There is good reason to believe that 
free expansion, while giving a good overall description for the data, is not a perfect model. 
Individual knots of ejecta become excited, and therefore visible, as they encounter the reverse 
shock, a process that also leads to the deceleration of these knots. Furthermore, some new 
material may become excited while that which has been shocked earlier fades, a process 
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that will also mimic deceleration. And finally, although ejecta knots are themselves moving 
radially, the reverse shock may have developed a non-spherical geometry in some directions, 
so knot-shock encounters can result in apparent changes in direction as well as decelerations 
for the excited filaments. All of these effects are seen in G292's "young er cousin" Gas A 



where the velocities are up to 5 times faste r, and the age 10 times younger (IFesen et al.ll2001 



Thorstensen et al.ll200ll : iFesen et al.ll2006l ). 



While it remains clear that the constant-expansion model is a good overall description 
for the 0-rich filaments in G292, the large values mean that we cannot reliably use 
incremental contours in space to determine confidence limit s on the model par ameters. As 



an alternative we turn to the bootstrap technique described in lPress et al.l (120071 ). Briefly, in 
this technique one starts with the original sample of M points (in our case, the 67 filaments), 
and randomly selects exactly M values, with replacement. Thus some points may be selected 
multiple times; others not at all. We fit the new data sample with the same 3-parameter 
model (Eq. 4) as the real data to determine {xo,yo) and r. We then repeat this process for 
100,000 trials, and consider the resulting ensemble of model parameters. The mean values 
from the 100,000 trials are in extremely close agreement with those from the unique best fit 
to the actual data, but now we can examine the frequency to give a robust estimate for the 
confidence limits. 

The resulting allowed region for the expansion center is shown in Fig. 7, an enlarge- 
ment of the central region of G292. The 90%-confidence contours measure approximately 
A(R.A.) =+J-2'', A(Dec.) = ±4.8". The uncertainty in age is ±60 yr (Id), or tloo yr (90%- 
confidence). 

We have investigated the differences between our measured values for the proper motions 
and the best-fit free-expansion model to check for a possible systematic pattern. In Fig. 8 we 
show the measured values as red vectors (identical to those in Fig. 4), and the corresponding 
proper motions from the model in blue. The uncertainties are indicated as ellipses at the 
ends of the measured-value vectors. While there are a number of cases where the deviation 
between measured and model is large compared with the uncertainty (thus producing the 
high value), we can discern no pattern to these. In several of these cases, the error ellipse 
is highly elongated, i.e., the formal uncertainty is much higher in one coordinate than the 
other. In at least some of them, the small uncertainty in one direction may have resulted 
from a fortuitously low dispersion among a small number of measurements. 

We have also carried out fits using the average proper-motion values calculated using 
the simpler linear-fit method of §3.3, and found very similar results: an expansion center 
within 2" of that given in Eq. 5, with somewhat larger but similarly shaped error contours. 
The age resulting from these fits was 2900 ± 80 yr (Icr), overlapping with the value from 
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Eq. 5 at the Icr level. Given the great similarity in the two sets of results, we shall use only 
those from our preferred multi-epoch averages (Eq. 5) in the remainder of this paper. 



Discussion 



G292 shows a somewhat irregular morphology at all wavelengths. The largest feature, 
and the one which comes closest to circular symmetry, is the outer "plateau" region, as 
shown in the 20-cn i radio image taken from the Australia Telescope Compact Array by 
Gaensler fc Wallace! (120031 ). This shows a sharp outer shell with radius ~ 4', centered at 
R.A. (2000.0) = ll'^24™34.8^ Dec. (2000.0) = -59°15'52.9", with an uncertainty of ±5" in 
both directions. The same image shows that the brightest radio emission is concentrated 
within a much brighter inner core of radius ~ 2^ with a cen ter displaced ~ 30" east and 



south from that of the outer shell. I Gaensler fc Wallace! (12003! ) identify the core as a pulsar- 
wind nebula associated with the PSR J1124-5916 (see below), and the plateau boundary as 
the outer shell where the supernova blast wave is interacting with its environment. 

Our measured expansion center (Eq. 5) is remarkably close to the center of this outer 
radio shell: only 3" away, as shown in Fig. 7. While such close agreement (within the 
uncertainties of both center determinations) may be fortuitous, this common center makes it 
appear likely that the outer blast wave from the explosion that has shaped the radio shell has 
approximate spherical symmetry. Certainly it is more symmetrical than the distribution of 
inner ejecta that we are now seeing as [O IIl]-emitting filaments (see subsequent discussion), 
or than the bright core of radio emission, or than the X-ray-emitting plasma , which conaprises 



mate rial with varying degrees of enrichment from the supernova ejecta (jPark et al. 
2007h . 



2004 



The young pulsar PSR J1124-5916, located near the center of G 292, is almost certainly 
the compact remnant from the same supernova (jCamilo et al.l 12002! ). This has already been 
underscored by the presence of a prominent pulsar-wind nebula in X -rays, discove r ed by 
Hughes et al.! (l200ll ) prior to the discovery of the relatively faint pulsar. ICamilo et al.l (12002! ) 
have measured the spin-down age of the pulsar at 2900 yr, quite close to the 2990-year 
expansion age for the G292 remnant and further cementing their association. 

The pulsar is now located at R.A. (2000.0) = ll'^24"39.1^ Dec. (2000.0) = -59°16'20", 
46" southeast of the expansion center (Eq. 5). If we assume that the pulsar and the oxygen- 
rich filaments were launched at the same time, about 3000 years ago, we would expect 
the pulsar to be moving southeastward at 16 masyr^^, or a transverse space velocity of 
440 (ig kms~^. While fast to be sure, this is well within the normal range for pulsar birth 
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velocities (ICaraved Il993l : iFrail et al.lll994l ). and far short of the transverse velocity of the 
unpulsed neutron star RX J0822-4300 inside the oxygen-rich SNR Pu ppis A, directly mea - 
sured (using X-ray data from Cha ndra) to be at least 10 00 kms~^ (IHui &: Beckerl |2006| ). 
and probably closer to 1600 kms~^ (jWinkler fc Petrd 120071 ). A similar direct proper-motion 
measurement for the pulsar in G292 will present an extreme challenge. Nevertheless, the 
displacement of the pulsar from the expansion center of the remnant is further evidence that 
core-collapse supernovae have significant asymmetries. 



Ghavamian et al.l (120051) used imaging Fabry-Perot spectroscopy to explore the kine- 
matics of G292 through measurement of radial velocities for the [O III]-emitting filaments 
in the central region of G292 — a field including all but the filaments most distant from the 
center (and thus with the highest transverse velocity). They found at least some filaments 
over the entire velocity range they investigated: —1440 to +1700 kms~^. They plotted the 
radial velocity of the filaments as a function of the ir projected distance from the geometric 
center of the remnant, measured by them from the iGaensler fc Wallace! (120031 ) radio images 
as a point about 22" north of the center we determined. These data show that the ejecta 
filaments are irregularly distributed around a thick ellipsoidal shell in the position-radial 
velocity plane, suggesting a shell expanding with a speed of ~ 1700 kms~^. A simple scaling 
argument, assuming that the shell is expanding with transverse expansion velocity equal to 
their observed radial velocity extremes, led them to an age estimate of (3000 - 3400) yr. 

Our own measurement of the age, based as it is on actual measurement of the transverse 
motions of the filaments rather than estimating them from a model, is both more reliable 
and more precise , and is also independent of the distance to G292. The fact that the 
Ghavamian et al.l ( 20051) age measurernent is roughly consistent with ours indicates that the 
distance of 6 kpc (IGaensler fc Wallace! l2003l ) is probably about right. 



Taken together, the proper-motion and radial-velocity measurements indicate a definite 
asymmetry in the distribution of ejecta-dominated filaments that are now visible in G292. 
The proper-motion measurements, and also the spatial distribution of oxygen-rich filaments, 
show that the ejecta along an axis oriented roughly N-S are moving significantly faster than 
those along the E-W axis. In the E-W direction, we find —63 < fi^ < 52 masyr^^, or 
— 1800 dQ < Vx < 1490 d^ kms"^, compared with —125 < < 82 mas yr~\ or — 3570 rfg < 
Vy < 2340 (ig kms~^ in the N-S direction. Furthermore. !Ghavamian et al.! (12005! ) measured 
radial velocities —1440 kms~^ < Vrad ^ 1700 kms~^, virtually identical with the E-W 
range we have measured, so it is clear that the fastest moving visible filaments are those 
moving north and southo The most straightforward interpretation of these results is that the 



iGhavamian et al.l ( 20051 ) observed some filaments throughout the radial velocity range they examined, 
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current distribution reflects an asymmetry in the supernova itself, and that core material was 
ejected fastest along an axis oriented approximately N-S in the plane of the sky. However, 
the alternative of circumstellar material, distributed roughly in a plane perpendicular to 
the N-S axis and shaping the flow of ejecta, cannot be excluded. An ejecta distribution 
similar to that in G292 is found in more extreme form in Gas A, where the fastest knots are 



located in a prominent jet to the NE , with a less obvious jet to the SW (IFesen et al.lll988 



Fesen fc Gundersonlll996l : lFesenll200ll ). Asymmetric, possibly bipolar explosions seem likely 
to have produced both of these oxygen-rich remnants, and perhaps this is a general feature 
of core-collapse supernovae. 

With proper-motion measurements for the oxygen-rich filaments in hand, it will be pos- 
sible to combine these with radial velocity measurements to determine a fuller 3-dimensional 
model for the distribution of the visible ejecta in G292. We have recently obtained spectra 
of dozens of filaments that can be used in conjunction with Fabry-Perot images to pursue 
this project, but this analysis is left until a subsequent paper. 



We gratefully acknowledge the outstanding support, typical of the mountain staff at 
GTIO, during the observations from numerous observing runs over two decades that yielded 
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with. This work has been supported financially by the NSF, through grant AST-0307613 
to P.F.W., and from NASA, through grant NAG 5-8020 to P.F.W. and Chandra grants 
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Table 1. Journal of CTIO Observations 









Scale 






Filter 








Date (U.T.) 


Telescope 


CCD 


(" pixel- 1) 


Field (') 


Line 


Ac (A) 


AA* (A) 


Exposure (s) 


Observers 


1986 Nov 29 


4.0 m 


RCA 


0.589 


5.0x2.9 


[O III] 


5020 


54 


2 X 1000 


PFW, R.P. 
Kirshner, 
J.R Hughes 


1991 Apr 19 


4.0 m 


TEK IK #1 


0.470 


8.0 


[O III] 
Blue 


5020 
4770 


54 
100 


5 X 600 
4x600 


PFW, KSL 


1999 Jan 20-23 


0.9 m 


TEK 2K #3 


0.401 


13.7 


[O III] 
Green 


5026 
5135 


39 
90 


3 X 1200 
3 X 600 


PFW, D. Paul 


2000 Jan 30-31 


0.9 m 


TEK2K #3 


0.401 


13.7 


[O III] 
Green 


5026 
5135 


39 
90 


2 X 1200 
1 X 1000 

2 X 600 
1 X 500 


PFW, KSL, 
E. Galle 


2002 Mar 20-23 


0.9 m 


TEK 2K #3 


0.401 


13.7 


[O III] 
Green 


5006 
5135 


60 
90 


5 X 1000 
5 X 500 


PFW, KSL, 
C. Reith 


2006 Mar 30, 
Apr 2 


0.9 m 


TEK 2K #3 


0.401 


13.7 


[0 III] 
Green 


5006 
5135 


60 
90 


4 X 1000 

3 X 800 

4 X 500 
3 X 400 


PFW, KSL, 
K. Twelker 


2008 Mar 4-6 


0.9 m 


TEK 2K #3 


0.401 


13.7 


[O III] 
Green 


5006 
5135 


60 
90 


5 X 1000 
5 X 500 


PFW, KSL 



^FuU width at half maximum in the telescope beam. 



Table 2. Proper Motions of Individual Filaments in the SNR G292.0+1.8 



Fil. No. R.A.^ Dec.'' (is cr^^ No. Epochs 

" " mas yr~^ mas yr~^ mas yr~^ mas yr"-*^ 



1 


-137.7 


-112.9 


-64.9 


16.6 


-18.5 


3.4 


6 


2 


-90.2 


6.4 


-60.0 


7.4 


12.6 


3.5 


4 


3 


-90.0 


-157.4 


-52.0 


4.9 


-54.7 


4.9 


5 


4 


-89.4 


53.5 


-34.4 


10.0 


24.7 


5.7 


5 


5 


-85.6 


-130.8 


-31.8 


7.0 


-49.5 


7.3 


6 


6 


-83.8 


-185.9 


-32.1 


5.6 


-57.7 


3.1 


6 


7 


-83.0 


-193.9 


-32.8 


3.8 


-58.1 


2.6 


5 


o 


-82.2 


-291.8 


-24.1 


3.9 


-79.9 


9.3 


b 


9 


-65.5 


-290.6 


-42.2 


4.9 


-122.1 


8.1 


5 


10 


-64.8 


-300.2 


-44.0 


5.6 


-105.5 


7.6 


5 


11 


-63.1 


-17.9 


-26.9 


8.0 


-19.2 


10.1 


4 


12 


-53.5 


-126.4 


-24.0 


3.2 


-26.5 


9.0 


6 


13 


-51.3 


-196.3 


-20.9 


4.2 


-61.2 


8.4 


6 


14 


-48.9 


125.3 


-19.7 


2.6 


46.7 


2.9 


6 


15 


-32.8 


98.4 


-4.2 


3.9 


42.8 


3.7 


6 


16 


-32.2 


22.8 


-20.0 


0.5 


9.5 


1.6 


4 


17 


-28.6 


-135.6 


-14.5 


7.6 


-22.5 


4.3 


5 


18 


-28.4 


-305.4 


-14.9 


3.9 


-101.2 


8.8 


5 


1 Q 




lO.O 


-J-tJ.U 






9 A 


A 


20 


-15.3 


-304.8 


-25.0 


6.6 


-107.1 


8.3 


5 


21 


-11.6 


-120.3 


-18.1 


6.1 


-19.4 


7.6 


4 


22 


-7.4 


91.6 


-13.6 


3.6 


35.1 


3.3 


4 


23 


6.7 


52.2 


-8.1 


0.7 


19.7 


3.7 


5 


24 


7.7 


113.7 


-17.0 


2.1 


23.8 


6.0 


6 


25 


7.7 


175.4 


2.3 


1.9 


62.4 


2.8 


4 


26 


9.7 


72.3 


-8.1 


2.0 


35.0 


3.5 


6 


27 


9.9 


121.1 


-7.8 


3.1 


40.8 


4.5 


4 


28 


15.8 


109.8 


2.6 


1.9 


39.2 


2.0 


4 


29 


19.7 


86.8 


3.2 


3.9 


32.2 


5.7 


6 


30 


35.2 


62.3 


12.6 


5.7 


32.2 


2.6 


5 


31 


36.9 


140.5 


5.6 


1.0 


54.8 


4.8 


4 



Table 2 — Continued 



Fil. No. R.A.*^ Dec.'^ //q, a^^ ns cr^^ No. Epochs 





II 


II 


mas yr ^ 


mas yr ^ 


mas yr ^ 


mas yr ^ 




32 


55.4 


171.2 


23.3 


2.4 


56.4 


4.7 


6 


33 


55.8 


194.8 


12.9 


2.7 


81.0 


2.5 


6 


34 


67.4 


174.6 


7.1 


3.3 


64.0 


8.0 


4 


35 


67.8 


-251.7 


33.4 


13.6 


-91.1 


5.7 


3 


36 


73.4 


-260.9 


9.1 


4.1 


-78.6 


6.8 


5 


37 


78.6 


-226.0 


27.4 


3.1 


-82.5 


6.5 


3 


38 


83.6 


-200.1 


20.3 


2.2 


-61.8 


2.5 


6 


39 


89.5 


-45.1 


20.6 


1.0 


-4.2 


1.3 


7 


40 


90.1 


-31.7 


15.4 


2.0 


-11.8 


2.3 


7 


41 


93.8 


-121.6 


20.5 


4.2 


-25.2 


3.9 


6 


42 


94.5 


-60.8 


26.4 


2.3 


-16.9 


1.1 


7 


43 


96.3 


-21.2 


21.1 


1.5 


6.1 


1.2 


7 


44 


97.6 


146.3 


29.0 


5.9 


67.0 


8.2 


4 


45 


104.7 


-6.3 


25.4 


1.7 


1.9 


1.5 


7 


46 


105.7 


-31.3 


22.5 


1.9 


-0.8 


1.6 


7 


47 


106.3 


-37.0 


26.5 


1.3 


-7.0 


1.4 


7 


48 


106.5 


-52.4 


20.1 


2.5 


-10.2 


0.8 


7 


49 


107.5 


-78.6 


24.1 


1.8 


-16.5 


1.9 


7 


50 


112.0 


-28.0 


27.4 


1.1 


-8.6 


1.9 


7 


51 


113.7 


-108.1 


30.9 


5.6 


-27.6 


3.0 


6 


52 


118.7 


-14.5 


40.3 


2.4 


10.0 


1.4 


7 


53 


122.3 


-97.5 


34.5 


2.6 


-25.9 


2.0 


7 


54 


123.7 


-33.5 


32.6 


3.2 


2.7 


1.9 


7 


55 


125.4 


-7.0 


33.9 


1.4 


0.2 


1.5 


7 


56 


127.4 


6.0 


36.7 


1.7 


11.6 


2.0 


7 


57 


131.8 


12.2 


32.5 


1.3 


19.3 


2.0 


7 


58 


132.5 


-18.1 


40.8 


1.9 


-0.6 


2.1 


7 


59 


135.0 


-50.8 


35.0 


2.6 


-10.2 


1.5 


7 


60 


135.8 


-11.3 


37.9 


1.0 


3.0 


1.5 


7 


61 


137.0 


-46.0 


43.6 


2.5 


-15.4 


1.8 


7 


62 


139.4 


63.1 


26.2 


1.9 


24.5 


0.9 


7 
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Table 2 — Continued 



Fil. No. 


R.A.'^ 


Dec* 






^J'5 




No. Epochs 




11 


II 


mas yr~^ 


mas yr~^ 


mas yr~^ 


mas yr~^ 




63 


142.2 


73.6 


33.5 


9.8 


35.2 


6.9 


6 


64 


144.3 


7.2 


50.0 


3.3 


16.8 


2.7 


7 


65 


153.0 


-13.2 


48.3 


2.5 


4.4 


0.7 


7 


66 


159.9 


48.8 


48.6 


3.8 


18.5 


2.1 


5 


67 


161.0 


58.5 


51.1 


3.3 


27.6 


1.8 


7 



^Offsets in arcsec relative to R.A.(J2000.) = 11 24 31.0, Dec.(J2000.) = - 59 15 30 (an 
arbitrarily chosen point near the remnant canter). Positive offsets are east and north. 
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Fig. 1.— Continuum-subtracted image of G292.0-hl.8 in the light of [O III] A 5007. The 
small blue boxes indicate the locations of 67 filaments whose proper motions have been 
measured. The larger boxes show the regions that appear in the enlarged difference images 
in Fig. 2. In order to show the fainter filaments clearly, many of the bright filaments in the 
eastern "spur" region (within and near box a) are saturated in this display, and may appear 
to blur together. In fact, each of the filaments we have used for our measurements (those in 
the small blue boxes) are distinct. 
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Fig. 2. — The top panels (a - c) are enlargements of the three I'-square sections indicated 
in Fig. 1, with one particularly prominent filament in each indicated by the arrow. We use 
these same filaments to illustrate our analysis methods in subsequent figures. The lower 
panels show the same sections as the difference between images taken at different epochs. 
In (d) the difference is between 1986 and 2008 images, while (e) and (f) use the difference 
between 1991 and 2008 images. (We have used the longest baseline available showing these 
filaments.) In each case emission at the earlier epoch appears as white, while the later epoch 
appears as black. Outward motion, while subtle, is apparent in each case. 
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Fig. 3. — Proper motion velocities (in x and y) of three prominent filaments. (These are 
the same filaments identified in Fig. 2.) Each epoch pair (21 pairs for filament 60, 15 for 
filaments 7 and 33) gives an independent measure of the motion. The dispersion is smallest 
for the measurements with the longest baselines, and naturally we weight these more heavily 
in determining the average of all measurements (solid lines) which we take as the best overall 
statistic (§3.2). 
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Fig. 4. — Continuum-subtracted image of G292, with vectors representing the measured 
proper motions of oxygen-rich filaments projected forward 1000 years at the rates we have 
determined. The vectors clearly radiate outward from a point near the center of the remnant, 
indicated by the cross (see §4; the central region is enlarged in Fig. 7). 
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Fig. 5. — Displacements (in x and y) of three prominent filaments (the same filaments as in 
Figs. 2 and 3) relative to their position at our 2002 epoch. The slope of the best linear fit 
gives a measure of the proper motion, as described in §3.3 
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Fig. 6. — Plots (in x and y) of the measured proper motions for 67 filaments in G292, as a 
function of position. (The position is measured in arcsec relative to an arbitrary zero point 
near the center of G292.) The error bars represent our best estimate of the uncertainties in 
measurements for individual filaments, according to the multiple-epoch technique described 
in §3.2. The data are in good overall agreement with a constant-velocity expansion model, 
represented by the best-fit lines. The inverse slope of the line, constrained to be the same 
for both the x- and y-component fits, gives the age, and the intercepts with fi = give the 
coordinates of the expansion center. 
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Fig. 7. — In this 2' square central section of an unsubtracted [O III] A 5007 image of G292, 
our expansion center is indicated by the red cross, with the 90%-confidence region also 
indicated in red. The blue cross indicates the position of PSR J1124-5916, 46" SE of th e 
expansion center. The center of the outer radio shell, as given by lGaensler fc Wallace! (12003! ) . 
indicated by the magenta cross, is only 3" from our expansion center. Some of the brighter 
[O III] filaments can be seen at the left (east) side of the figure. 
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Fig. 8. — The measured proper motions projected ahead 1000 yr are again represented as red 
vectors on the image of G292, exactly as in Fig. 4, and in addition the 1000- yr proper motions 
for the best-fit model are represented as blue vectors. Uncertainties in the measured values 
are represented by ellipses at the ends of the measured-value vectors. The cross represents 
the best-fit expansion center, identical to that in Figs. 4 and 7. While many of the model 
vectors lie well outside the corresponding error ellipse, there seems to be no strong pattern 
to the deviations. 



